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Abstract 

We explore the anomaly induced effective QCD meson potential in the 
framework of the effective Lagrangian approach. We suggest a decoupling 
procedure, when a flavored quark becomes massive, which mimics the one 
employed by Seiberg for supersymmetric gauge theories. It is seen that, after 
decoupling, the QCD potential naturally converts to the one with one less 
flavor. We study the N c and Nf dependence of the 7/ mass. 
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I. INTRODUCTION 



In the last few years there has been a great flurry of interest in the effective Lagrangian 
approach to supersymmetric gauge theories. This was stimulated by some papers of Seiberg 
|]J and Seiberg and Witten in which a number of fascinating "exact results" were ob- 
tained. In particular Seiberg has provided a very interesting picture for different phases of 
supersymmetric gauge theories. There are already several interesting review articles 

It is natural to hope that information obtained from the more highly constrained super- 
symmetric gauge theories can be used to learn more about ordinary gauge theories, notably 
QCD. In a recent paper || it was shown how the effective Lagrangian for super QCD might 
go over to the one for ordinary QCD by a mechanism whereby the gluinos and squarks in 
the fundamental theory decouple below a given supersymmetric breaking scale m. To ac- 
complish this goal a suitable choice of possible supersymmetry breaking terms was used. An 
amusing feature of the model was the emergence of the ordinary QCD degrees of freedom 
which had been hidden in the auxiliary fields of the supersymmetric effective Lagrangian. 
Constraints on the supersymmetry breaking terms were obtained by requiring the trace of 
the energy momentum tensor to agree (at one loop level), once supersymmetry was broken, 
with that of ordinary QCD. It was also noticed that a reasonable initial picture resulted 
by neglecting the Kahler terms in the original supersymmetric effective Lagrangian. This 
feature is analogous to Seiberg's treatment [I]] of supersymmetric effective Lagrangians with 
different numbers of flavors. It led to a dominant piece of the QCD effective Lagrangian 
possessing a kind of tree level "holomorphicity" . Physically this corresponds to the explicit 
realization of the axial and trace anomalies by the model. 

In this letter we will further explore the " holomorphic" part of the potential for QCD 
with Nf(< N c ) massless quarks as obtained by breaking super QCD according to the scheme 
above (see section III of 0). In particular, we shall study the analog of the decoupling 
procedure used in |TJ when one flavor becomes heavy. In the present case we no longer have 
the protection of supersymmetry for deriving exact results so the procedure will be carried 
out at the one loop level. The analog of the Affleck-Dine-Seiberg (ADS) superpotential @ 
turns out to be |1 the holomorphic part of the potential: 



V(M) = —C (N c , Nf) 



detM 



ll(N c -Nf) 

+ h.c. , (1.1) 



where A is the invariant QCD scale and the meson matrix M contains Nj scalars and Nj 
pseudoscalars. A particular choice of the dimensionless positive quantity C (N c , Nf) was 



1 



made in |J. Here we find the constraints on C (JV C , Nf) which follow from decoupling. 
The result can be used to suggest that the well known []8||| large N c behavior of the 
(pseudoscalar singlet) meson mass should also include an Nf dependence of the form: 

M v' K F^V A2 {Nf K Nc) ■ (L2) 

This could be considered an indication that, while M v i is suppressed in the large N c limit, 
there is an enhancement mechanism for Nf near N c — \ to explain its relatively large nu- 
merical value. 



II. MESON POTENTIAL AND QUARK DECOUPLING 

The various super QCD (SQCD) effective Lagrangians are essentially derived using 
the anomaly structure of the underlying gauge theory as well as supersymmetry. It is thus 
interesting to first see how much of the effective QCD Lagrangians we get after decoupling 
the super-partner fields, follows just from anomaly considerations. In ordinary QCD the 
two relevant anomalies are the trace and Adler-Bell-Jackiw Ua{1) anomalies: 

02 = -b^- 2 FrF mn ., a = 2bH , (2.1) 
Su a (i)Vqcd = N f a^- 2 e mnrs FrF r a s = AN f aG . (2.2) 
11 2 

Nf is the number of flavors, b = — N c Nf is the coefficient of the one loop beta function 

3 3 

and left handed quarks have unit axial charge. H and G can be interpreted as composite 
operators describing, in the confining regime, the scalar and pseudoscalar glueball fields |TU 



The general effective (i.e. in terms of composite operators) potential which, at tree level, 
saturates the one loop anomalies is 



V = -bFj2 C ^n(^)+h.c.+V I , (2.3) 



n 



where O n is an operator built out of the relevant degrees of freedom at a given scale \x 
(glueballs, mesons, baryons, ..) with mass dimension n and axial charge q n , F = H + i5G 
and Vi is a scale and £7a(1) invariant potential. The anomaly constraints are: 

5>=1, E%„ = ^. (2.4) 

A is connected with the invariant scale of the theory and the 9 parameter (at one loop) via 
A b = [i b e~ T with r = -^^y — id. At high scales ( ^> A) we expect to approach the classical 
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regime where the glue potential displays the holomorphic structure —r(H + iG) +h.c This 
suggests choosing 5 = 1 and hence F = H + %G. Here we assume that the lightest meson 
(M)R and glueball (H and G) fields are the relevant degrees of freedom. It is amusing to 
notice that the field associated with the Ua(\) transformation (77') can only appear in the 
first term of the potential in Eq. ( |2.3| ). 

We can arrange the anomaly potential (first term in Eq. ( |2.3| )) to be "holomorphic" by 
setting O n = O n (F, detM). By holomorphic we mean a potential of the form + h.c. 
where x is a function of the generic complex field $. Now, in Ref. || it was shown that a 
suitable decoupling of the holomorphic superpotential in the Taylor Veneziano Yankielowicz 
(TVY) model [[?[ yielded an important "holomorphic" contribution to the ordinary potential 
for QCD. For Nf(< N c ) flavors the result was of the form 



V(F,M) = -^(N c -N f )F 



In 



-N C F" 
7A 4 , 



Fin 



/ detM \ 



+A (N f ,N c )F + h.c. 



(2.5) 



where 7 



\2N r . - AN 



f 



11N - 2N 



f 



and A (N c , Nf) is a dimensionless constant which cannot be fixed by 



saturating the QCD anomalies and assuming holomorphicity. This is because the term AF 



does not contribute to 9™ and is also a chiral singlet. The potential in Eq. (2J3) is easily 

seen to be consistent with the general form in Eq. (|2.3|) , and the associated constraints in 

Eq. ( p.4| ) when the scale dimension of M is fixed to be 3 and 5 = 1. (The constants are 
N f N f N 

c 4 



N f , Nf 
1 — 3—— and c^N r — 



b J1,/ b 

By integrating out the "heavy" glueball field F one obtains the "light" meson field 
potential for M. This is given in Eq. (|1.1| ) wherein the unknown coefficient C (N c , Nf) is 
related to A(N C , Nf) by the following expression: 



A(N c ,N f ) 



11 
12 



(N c -N f )ln 



12 



Nr. 



11 1 (N c -N f ) 



C{N c ,N f ) 



(2.6) 



Similar effective potential models and some related phenomenological questions have already 



been discussed in the literature IJTTHI5 



The potential for the meson variables in Eq. ( |1 . 1| ) is similar to the effective ADS super- 
potential for massless super QCD theory with Nf < N c 0] 



W 



ADS 



(T) = - (JV C - N f ) 



A 



3N C -N f 



detT 



(2.7) 



*In Ref. §§ M ~ —qnqi was denoted Ft- 
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where T ~ QQ is the composite meson superfield, Q and Q are the quark super fields and As 

is the invariant scale of SQCD. In the instanton generated super potential the exponent of 

As is the coefficient of the super symmetric beta function. In the potential for the ordinary 

meson variables Eq. ( |1 . 1|) the exponent inside the square brackets of the QCD invariant scale 

11 2 

A is similarly the coefficient of the one loop QCD beta function b = — iV c Nf. As in 

the SQCD case the ordinary QCD potential Eq. ( |1.1| ) can also be constructed if we assign 
to A b charge 2Nf under the matter field axial transformation, as prescribed by the #-angle 
variation, and if holomorphicity in the coupling constant is assumed. 

An important difference with respect to the ADS superpotential, as already explained 
in Ref. pi, is the fact that the QCD potential displays a fall to the origin while the ADS 
potential does not have a minimum for finite values of the squark condensate. In || it was 
shown that the fall to the origin can be cured by introducing non holomorphic terms which 
also trigger spontaneous chiral symmetry breaking. 

The mesonic holomorphic potential, due to the presence of the appropriate beta function, 
is expected to hold for any Nf < N c . In particular a non trivial check would be to decouple 
a single heavy flavor and show that the potential reduces to the potential for one less flavor 
fermion. A similar procedure, at the supersymmetric level, was used by Seiberg to check 



the validity of the ADS superpotential. This is in the same spirit as the well known fl6 
criterion for decoupling a heavy flavor (at the one loop level). For a small mass m of the Nf 
-th quark the perturbative contribution, prescribed in the fundamental lagrangian, of the 
mass operator to the potential is 



V(M) = ruM^l + h.c. . (2.8) 

To achieve a complete decoupling we generalize the perturbative mass operator to 

V m = -m A Mxf + h.c. , (2.9) 

where dimensional analysis requires A = 4 — 3D We interpret the departure from unity of 
A as an effective dynamical evolution of the chiral symmetry breaking operator for large 
values of m. The total potential for large m, 



V T {M) = -C(N c ,N f ) 



AT N c-l N f 
detM 



if N c -Nf) , I ■ 



— m A Mfr +h.c. , (2.10) 



where A is the appropriate invariant scale for Nf flavors and N c colors, must convert to the 
potential for Nf — 1 light quarks. Notice that the generalized mass term preserves the holo- 
morphic structure. Here it is being assumed that the non-holomorphic terms in the potential 
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can be treated as small pertubations. A similar mass operator was already introduced in 
Ref. H to appropriately decouple the gluino in the underlying theory. Eliminating the heavy 
degrees of freedom by their equations of motion and substituting back in the potential gives 



where w 



V(M) = - 
12 



T + w 



w 



C (N c , N f ) w 



r 

r+w 



A b(N c ,N f ) m i 

detM 



r+w 



+ h.c. , 



(2.11) 



, , b(N c , Nf) is the coefficient of the beta function for the theory with 
11 (N c -Nf) K SJ 

N c colors and Nf flavors and M is the meson matrix for Nf — 1 flavors. In the standard 
physical picture the gauge coupling constant evolves according to the QCD beta-function 
for Nf flavors above scale m and according to the QCD beta-function for Nf — 1 flavors 

.. e ; 

b(N c ,N f ) 



57T 



exp 



the matching at jj, — m requires 



'A 



N c ,N f 

m 



\ \ b(N c ,N f -l) 

m j 



which yields 



A 



b(Nc,N f -l) _ A b(N c ,N f )^2 



N c ,N f -l 



N C ,N 



f 



(2.12) 



If we require the presence of the chiral symmetry breaking term in Eq. ( |2.10| ) to convert 

2 

the Nf theory into the Nf — 1 theory, the ratio A/T is fixed to - and by using the relation 

3 

A = 4 — 3r one derives 
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A = — , 
11' 



11 



(2.13) 



(2.14) 



This valueQ of F is exactly what is needed to get the correct exponent 

Tw _ 12 
T + w ~ 11 (N c -N f + 1) ' 

in Eq. ( |2.11| ). A consistent decoupling is obtained provided that the coefficient C (N c ,Nf) 
obeys the following recursive relation 



\C(N c ,N f y 


N c -Nf 


N c — Nf 





C(N c ,N f -l) 
N c — Nf + 1 



Nc-N f +1 



(2.15) 



The general solution of Eq. (|2.15 ) is finally. 



t These are the same exponents found in section II of || for decoupling a massive gluino. Here 
we can also use, as in ||, the one loop matching of the anomaly of the traced energy momentum 
tensor to derive the same exponents. 
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C(N c ,N f ) = (N c -N f )D(N c ) N <- N f . (2.16) 
Using Eq. ( |2.6|) , A (N c , Nf) can be expressed as the function of D: 

A (N c , N f ) = H (N c - N f ) In + ^\nD (N c ) . (2.17) 

It is interesting to contrast the result in Eq. ( |2.16|) for the coefficient of the "holomorphic" 
part of the QCD potential Eq. Q with Seiberg's result [Q C (N c , N f ) = (N c - Nf) for 
the coefficient of the ADS superpotential in Eq. ( [2.71 ). Clearly, in the SUSY case the analog 
of D(N C ) is just a constant (which turns out, in fact, to be unity by an instanton analysis 
0). This feature arises in the SUSY case because of the existence of squark fields which can 
break the gauge and flavor symmetries by the Higgs mechanism; it leads to C (N c , Nf) = 
C (N c — Nf). Since there are no appropriate scalar fields in ordinary QCD we do not expect 
such a feature. The possibility of a non-constant D(N C ) factor can thus be taken as an 
indication that there is no Higgs mechanism present. 



III. PHYSICAL EFFECTS 



As mentioned above, we must add to the "holomorphic" piece of the potential (Eq. ( |1 . 1| ) 
with Eq. ( |2.6| )) a " non- holomorphic" piece in order to prevent a "fall to the origin". This 
would also yield a non- zero value for (M„). The non-holomorphic piece will be required to 
neither contribute to the trace anomaly nor to the axial anomaly. At the level of Eq. fl2.3Q , 
where the "heavy" glueball fields have not yet been integrated out, such a potential could 
be of the form 



V I = ^ i A n (F*F)2~4™ Tr Um^M 



(3.1) 



where the A's are real constants. (For definiteness the terms have been restricted to those 
which are leading in the large N c limit). For our initial analysis we assume that Vj is small 
compared to the holomorphic piece and can be treated as a pertubation. An example of a 
single term with n > | was presented in Ref. 0. An analogous term without the heavy F 
fields is 



V/ = ]T4jr[(MtM) n 



(3.2) 



which should be added to Eq. (LI). The new feature of the present model is that the 
piece which mocks up the anomalies can be thought of as inheriting a specific holomorphic 
structure from the SQCD theory. 
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Since the model of Eq. ( |1.1|) plus Eq. (|3.2|) contains a great deal of arbitrariness it is 
natural to wonder what can be learned from it. The key point is that Eq. ( |3.2| ) possesses an 
axial U (Nf) L x U {Nf) R symmetry and cannot contribute to the rj' (pseudoscalar singlet) 
mass. Thus we may learn something about this quantity from Eq. ( |1 . 1|) directly. Now the 
rf field may be isolated from the meson matrix M as 



M = KU exp 



•q 



aAJNf 



(3.3) 



where K = K\ = U" 1 , detU = 1 and a is a dimensionless real constant. K contains the 
scalar fields and U the non flavor-singlet pseudoscalars. For our present purpose we may set 
K = (3A 3 1, where (3 is another dimensionless real constant. For simplicity we shall make 
the approximation a = (3 = 1, which will not affect our results in the large N c limit. Then 
substituting Eq. (|3.3| ) into Eq. fll.l]) , using Eq. ( j2.16Q and expanding the result to quadratic 
order in rj' yields the mass 



, 2N f /12\ 



2 



Mi = j^—^- {-) AD (N e )«°-«, . (3.4) 

Here it has been assumed that the rj' has the canonical Lagrangian mass term — -d m rj'd m rj'. 

If the factors a and j3 were to be included the right hand side of Eq. (|3.4|) would get multiplied 

12N f 

by oT 2 ft u(N c -N f ) ^ a ^ oeg no ^. jjg^g an y i ar g e jy c dependence and that of ft gets washed 
out. 

D(N C ) is unspecified by our model so we cannot predict M v i. We can however check the 

dependence on iV c of the quantities in Eq. (|3.4| ) for large N c . It is well known that 

Mj?/ ~ 1/N C for large A^ c . Furthermore, (F) (defined after Eq. ( |2.3| )) should behave like A^ c 

for large A^ c by standard counting arguments. Hence (see the discussion after Eq. (2.18) of 

i i 

Ref. II) A is expected to behave as (iV c (F)) 4 ~ A^ 2 at large N c . For consistency we thus 
expect the large A^ c behavior 

.Nr. 



D(N e )~[±-\ . (3.5) 



It is amusing to observe that when Nf is close to N c the resulting pole in Eq. ( |3.4| ) suggests 
a possible enhancement mechanism for the rj' mass. This could explain the unusually large 
value of this quantity in the realistic three flavor case. 

Of course the Nf = N c case is rather non-trivial. Equation ( |1.1| ) shows that the overall 
coefficient of our model, C (N c , Nf) should vanish for Nf = N c . This is what happens 
in the SQCD situation. There, Seiberg H shows that the effective superpotential should 
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be replaced by a quantum constraint which also involves the baryon superfields. Before 
considering the analog of these results in our context, let us try to understand the situation 
a little better by going back to the holomorphic part of the potential in Eq. ( |2.5| ), in which 
the heavy glueball field F has not yet been integrated out. Focusing on the pieces of Eq. ( p.5|) 
which are related to generating an rf mass we get in a schematic form 

V (F, M) ~ — (iV c - Nf) G 2 + 7]'G , (3.6) 



where G is the pseudoscalar glueball field defined in Eq. ( p2.2[ ) and inessential factors have 
been dropped. The first term is a wrong sign (for Nf < N c ) mass term for G and is obtained 
by expanding the first term of Eq. (|2.5| ). Integrating out G according to the equation of 
motion results in a correct sign mass term for the rf as long as Nf < N c . It is clear that 
this mechanism will not work for Nf > N c From our present treatment we can see what 
is likely to be happening in the case of ordinary QCD. Here the non-holomorphic piece in 
Eq. ( |3.1|) may contain a — G 2 piece when it is expanded and may thus generate an rf mass. 
We have been treating Vj in Eq. (|3.1|) as a small pertubation but, for consistency, it must be 
non- negligible especially for Nf > N c . The usual phenomenology is consistent with the large 
N c approximation to QCD so perhaps it is easiest to use the present model in that case too. 
In that event we would always have Nf <C N c . However it is certainly of interest to examine 
the small N c case. The present model suggests that something different is occurring as we 



let Nf vary near iV c . It is noteworthy that recent preliminary computer simulations [17 
indicate that there is a rather noticeable change in the quark condensate for Nf near N c . 

Now let us return to the discussion of the holomorphic part of the potential and point 
out that formally the Nf = N c case can be handled as in SQCD. When we put Nf = N c in 

Eq. rra 



V(F,M) - -Fin [ eMrf ^ A3N% ) + h.c. . (3.7) 



Integrating out F then gives the constraint on the field M, 

detM = e A ^' N ^A 3 ^ Nc , (3.8) 

rather than a potential V(M). In fact the holomorphic potential vanishes as we argued. 
Nevertheless this vanishing potential is consistent (as in the SQCD case) with decoupling 
to the Nf — N c — 1 theory. Following Eq. ( |2.9| ) we consider the N c = Nf potential with a 
heavy iVj-th quark of mass m: 

V = - m A [M^f + h.c. , (3.9) 
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where A = 8/11 and T = 12/11 as before. Now the constraint in Eq. ( |3.8j ) enables us to 
write 

M^fdetM - e A ^ N ^Af ^ N , (3.10) 

J C) c 

where M is the meson matrix for N c — 1 flavors. Substituting this back into Eq. (|3.9|) yields 

(A 3 jv \ r 
ffl y +/i.c. . (3.H) 
detM ) 

When one uses Eq. ( |2.12| ) and Eq. ( p. 17] ) this is recognized as precisely Eq. ( |LT| ) together 
with Eq. ( |2.16| ) for the Nf — N c — 1 case. This argument may also serve to clarify what is 
"happening" in the SQCD situation. 

A final question concerns the role of baryons when Nf = N c . Baryon superfields were 
shown to play an important role [|IJ for Nf = N c in SQCD. One might hope to decouple 
the superpartners of the baryon superfields B and B [|l]] to obtain information on the QCD 
three quark baryons. However it is easy to see that none of the components of B and B 
contain any three quark composite fields. Thus the quick answer would be that the baryon 
fields are irrelevant for the ordinary QCD effective Lagrangian. This point of view is perhaps 
reinforced by the usual treatment (for large iV c ) of baryons as solitons of the effective meson 
Lagrangian. Nevertheless, the arguments in the first part of section II show that ordinary 
baryon fields may readily be added to the holomorphic potential in a manner consistent with 
the anomaly structure of QCD. They may very well be relevant for a small N c treatment of 
the theory. We hope to return to this question elsewhere. 
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